Coexistence of a nearly spin-singlet state and antiferromagnetic long-range order in 

quantum spin system Cu 2 CdB 2 06 
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A 1/2 magnetization plateau in magnetic fields above 23 T and antiferromagnetic (AF) long-range 
order (AFLRO) in low fields were found in Cu2CdB2 0e. Experimental results agree with quantum 
Monte Carlo results for an expected spin system. There are two kinds of Cu sites [Cu(l) and Cu(2)], 
which are located adjacent to each other. Unexpectedly, spins on the Cu(l) and Cu(2) sites are in 
a nearly spin-singlet state and form AFLRO, respectively, although interactions between the Cu(l) 
and Cu(2) spins cannot be ignored. Cu2CdB20e is the first material which shows such coexistence 
in an atomic scale. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.30.Cr 
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Coexistence or competition of plural different states 
has attracted much attention in condensed matter 
physics. Examples of such coexistence are an ordered 
stripe phase pj and spatial distribution of electronic en- 
ergy gaps Qj in high-T c cuprate superconductors, a mix- 
ture of insulating regions with charge ordering and metal- 
lic ferromagnetic domains in manganites and coexis- 
tence of superconductivity and antiferromagnetic order 
in i?Ni2B2C (R — Ho or Dy) These phenomena re- 
sult from strong correlations among electrons or spins, 
and therefore similar phenomena can be seen in various 
fields of condensed matter. 

Let us now consider whether coexistence of plural dif- 
ferent states occurs in quantum spin systems consisting 
of localized spins or not. In quantum spin systems, var- 
ious states are realized such as a spin-singlet state with 
a spin gap and a magnetically long-range ordered state. 
The former state is caused by quantum-mechanical ef- 
fects, while the latter can be usually understood in semi- 
classical pictures. These two states are representatives 
in quantum spin systems, and quite different from each 
other. There are examples of spin states, which possess 
characters of a spin-singlet state with a spin gap and 
antiferromagnetic (AF) long-range order (AFLRO). In 
spin-gap systems such as the spin-Peierls system CuGe03 
and the two- leg ladder system SrCu203 doped 
with impurities the ground state shows AFLRO, 

and two kinds of magnetic excitations exist. One is a 
high-energy excitation corresponding to a singlet-triple 
gap in pure systems, and the other is a low-energy spin- 
wave excitation reflecting the existence of AFLRO. How- 
ever, Cu sites of nearly spin-singlet states and those in 
AFLRO cannot be distinguished. Thus, the spin states 
in the spin-gap systems doped with impurities are differ- 
ent from the above-mentioned examples of coexistence. 
There is also an example of a magnet including both spin- 
singlet states and AFLRO which are almost decoupled. 
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In NH4CUCI3, the spin system consists of three differ- 
ent 5=1/2 dimer subsystems 0- States in two dimer 
subsystems are spin singlet, whereas AFLRO appears in 
another subsystem in low fields. Exchange interactions 
between spins in the spin-singlet states and those form- 
ing AFLRO are small in comparison with major interac- 
tions, meaning that the spin-singlet states and AFLRO 
are almost independent. Therefore, the coexistence of 
spin-singlet states and AFLRO in NH4CUCI3 can be un- 
derstood easily. 

At present, no magnet satisfying the following situa- 
tion, which cannot be easily expected, has been reported 
so far. Spins on some sites are in a nearly spin-singlet 
state and those on the other sites exhibit AFLRO, al- 
though exchange interactions between these two kinds of 
spins cannot be ignored. In such a case, it is natural to 
consider that either a spin-singlet state or AFLRO ap- 
pears [ic|. 

Unexpectedly, we found coexistence of a nearly spin- 
singlet state and AFLRO in C^CdB^Og. As is shown 
in Fig. 1, there are two kinds of Cu sites [Cu(l) and 
Cu(2)] in this cuprate In low magnetic fields, spins 
on the Cu(l) sites are in a nearly spin-singlet state, and 
those on the Cu(2) sites form AFLRO, although inter- 
actions between the Cu(l) and Cu(2) spins cannot be 
ignored. C^CdB^Oe is the first material showing such 
coexistence in quantum spin systems. 

We show an expected spin system of C^CdB^Og in 
Fig. 1. Taking the composition and reasonable valence 
of each ion (Cd 2+ , B 3+ , and 2 ~) into account, all Cu 
ions are divalent and have S = 1/2. Their positions are 
shown schematically in Fig. 1(a). Two crystallographic 
Cu sites [Cu(l) and Cu(2)] exist there along with five 
kinds of short Cu-Cu bonds, which are summarized in 
Table I. From Cu-O-Cu angles, signs of exchange inter- 
actions in Cu-Cu bonds 1, 2, and 3 are presumed to be 
AF. Values of the exchange interactions are defined as 
J\, J2, and J3, respectively. In bonds 4 and 5, no neigh- 
boring oxygen site exists to which both two Cu sites in 
each bond are connected. For that reason, exchange in- 
teractions in these bonds are expected to be smaller than 
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Ji, J 2 , and J3. In addition, Cu-Cu distances in the other 
bonds are larger than 4.64 A. Consequently, Cu 2 CdB 2 6 
is probably as the first approximation a compound in- 
cluding a unique spin system formed by J\, J 2 , and J3 
represented in Fig. 1(b). 

Crystalline powder of Cu2CdB 2 06 was synthesized 
using solid-state reaction method. We used X-ray 
diffraction measurement to confirm the formation of 
Cu2CdB206 and the absence of other materials. Us- 
ing an inductively coupled plasma atomic emission spec- 
trometer (ICP-AES), we estimated the composition of 
our sample as Cu1.9sCdB2.01O6.14, which closely resem- 
bled Cu 2 CdB20g. We measured magnetizations up to 
5 T using a superconducting quantum interference de- 
vice (SQUID) magnetometer produced by Quantum De- 
sign. High-field magnetizations up to 30 T were mea- 
sured using an extraction-type magnetometer in a hybrid 
magnet at the High Magnetic Field Center, NIMS. Spe- 
cific heat was measured using relaxation technique with 
PPMS, which was produced by Quantum Design. We 
calculated susceptibility and magnetization of the spin 
system in Fig. 1(b) by a quantum Monte Carlo (QMC) 
technique using the directed-loop algorithm in the path- 
integral formulation [l^ ■ The number of sites and Monte 
Carlo samples in the QMC simulations are one thousand 
and about one million, respectively. Finite-size effects 
and statistical errors are negligible in the scale of figures 
represented in this Letter. 

The solid curves in Figs. 2(a) and 2(b) represent tem- 
perature (T) dependence of magnetic susceptibility x{T) 
of Cu 2 CdB 2 6 in the magnetic field of H = 0.1 T. The 
susceptibility at high T obeys the Curie- Weiss law, which 
is expressed as C / (T + 0). Here C and 6 are a Curie con- 
stant and Weiss temperature, respectively. The inset of 
Fig. 2(a) shows inverse susceptibility. We roughly esti- 
mated C as 0.448 emu K/Cu mol from a slope of the in- 
verse susceptibility above 180 K and 6 as 47.4 K from the 
intersection of the dashed line and the T axis. We deter- 
mined independently the powder-averaged gyromagnetic 
ratio of Cu 2+ as g — 2.10 in an electron spin resonance 
(ESR) measurement at room temperature. Thereby, a 
value for C was calculated as C — 0.414 emu K/Cu mol 
assuming that all the Cu ions have S = 1/2. The two 
C values are similar to each other. This result indicates 
that almost all the Cu ions have S = 1/2. 

The susceptibility of Cu 2 CdB 2 0g has a maximum 
around 11 K and seems to reach a finite value at K. The 
inset of Fig. 2(b) represents dx(T)/dT. A lambda-type 
behavior is observed around 10 K. Figure 2(c) shows spe- 
cific heat C(T) in T, which has a peak at 9.8 K. Taking 
the expected signs of the exchange interactions into ac- 
count, AFLRO probably appears at Tn = 9.8 K. Besides, 
we found a spin-flop transition around 1.5 T in magne- 
tization M(H) measured by the SQUID magnetometer 
at 2.9 K. A small value of a spin-flop field -ffgp ~ 1.5 
T is reasonable because magnetic anisotropy of spins on 
Cu 2+ ions is usually small. For example, Hgp is about 
1 T in Cuo.96Zn .o4Ge0 3 [l|. The dotted curve in Fig. 



2(a) shows the QMC result of x(T) for the spin system 
in Fig. 1(b). It is explained later. 

The solid curve in Fig. 3 represents M(H) at 2.9 
K. The most prominent feature is a 1/2 magnetization 
plateau above 23 T. Diamonds represent the QMC re- 
sult of M(H) at 2.9 K for the spin system in Fig. 1(b) 
with Ji = 160, J 2 = 38.8, and J3 = 9.7 K and can 
reproduce well the experimental M(H). The QMC re- 
sult of x{T) for the same model is indicated by the 
dotted curve in Fig. 2(a) and agrees with the experi- 
mental x{T)- Therefore, the spin system in Fig. 1(b) 
can explain the magnetic properties of Cu 2 CdB 2 06 at 
least as the first approximation. As mentioned above, 
however, AFLRO appears at low T, which is caused by 
weak three-dimensional (3D) couplings ignored in the 
spin model. The discrepancy between experimental x(T) 
and the QMC result of x(^) a t low T is due to the ap- 
pearance of AFLRO. 

Spin states can be understood as follows (Fig. 4). The 
appearance of the 1 /2 magnetization plateau means that 
half of spins are in a nearly spin-singlet state with a dimer 
gap, which corresponds to the magnetic field to break the 
nearly spin-singlet state, and that the other spins are al- 
most polarized parallel to the applied field. Therefore, 
the magnetization cannot increase in the plateau region 
with an increase in the magnetic field. The QMC results 
in Fig. 3 show that moments of the Cu(l) spins (open 
circles) are very small in the low-field regime, indicating 
that the Cu(l) spins form nearly singlet pairs with the 
dimer gap below an end field of the plateau. The inset 
of Fig. 3 shows the QMC results of M(H) for a wider 
region of H. The magnetization starts to increase again 
around 108 T, which means that a value of the dimer gap 
is about 152 K for g = 2.1. The dimer gap corresponds to 
binding energy of singlet pairs and is mainly determined 
by the J\ interaction. Thus, the value of the dimer gap 
(152 K) is close to J x (160 K). The QMC results in Fig. 
3 show that moments of the Cu(2) spins (open squares) 
are almost polarized in the plateau region. Since the 
moments of the Cu(l) spins are small below 108 T, an 
effective interaction between two Cu(2) spins through the 
AF dimer formed by the Cu(l) spins is small. Thus, the 
Cu(2) spins can be almost saturated above 23 T in spite 
of the large values of J\ and J 2 [Fig. 4(b)]. AFLRO 
can appear in the Cu(2) spins [Fig. 4(a)] by this effec- 
tive interaction, the J3 interaction, and other weak 3D 
couplings which are omitted in our model |14|. The dis- 
tance of the nearest-neighbor Cu(l)-Cu(2) bond (bond 2) 
is only 3.22 A and the interaction in this bond (38.8 K) 
cannot be ignored in comparison with other interactions. 
Nevertheless, the coexistence of the nearly-singlet state 
of the Cu(l) spins and AFLRO of the Cu(2) spins ap- 
pears. We can say that this is "coexistence in an atomic 
scale" and such coexistence in quantum spin systems has 
not been found before the present study. 

As was described, both the Cu(l) and Cu(2) spins 
have non-zero moments. The average value of the to- 
tal spins is, on the contrary, exactly S = 0.25 in the 
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TABLE I: Interatomic distances and angles in Cu2CdB2 0e 



Cu-Cu distance (A) Cu-O-Cu path angle (°) Cu-0 distance (A) number of path sign 

bond 1 (Ji) 2.98 Cu(l)-Cu(l) Cu(l)-0(3)-Cu(l) 99^24 1.96, 1.96 2 AF~ 

bond 2 (J 2 ) 3.22 Cu(l)-Cu(2) Cu(l)-0(2)-Cu(2) 118.48 1.89, 1.86 1 AF 

bond 3 (J 3 ) 3.40 Cu(2)-Cu(2) Cu(2)-0(1)-Cu(2) 102.39 1.95, 2.40 1 AF 

bond 4 3.40 Cu(l)-Cu(l) none 

bond 5 3.56 Cu(l)-Cu(l) none 



plateau region This result indicates that a simple idea, 
in which some spins are perfectly spin singlet and oth- 
ers are fully polarized, cannot explain the magnetization 
plateau in C^CdP^Oe- Similar spin states must appear 
in other spin systems possessing magnetization plateaus. 
It is proved that magnetization curves at K may have 
plateaus at ml satisfying the formula of n(S' — m!) = in- 
teger 15] . Here, n is period of ground state, and S' and 
m' are the total spin and magnetization per unit cell, re- 
spectively. In our model, n = 1, S' = 2 (four 5=1/2 
spins), and the plateau appears at m' = 1. Thus, our re- 
sult is consistent with the theorem of Ref . 0] . We also 
calculated magnetization curves at 5.5 and 15 K and con- 
firmed that the values of magnetization at the plateau are 
also S = 0.25. 

In conclusion, we observed the 1/2 magnetization 
plateau in magnetic fields above 23 T and the antifer- 
romagnetic (AF) long-range order (AFLRO) in low mag- 
netic fields in Cu2CdB206- The experimental results ex- 
plained herein are consistent with the quantum Monte 
Carlo results for the expected spin system. The two kinds 
of Cu sites [Cu(l) and Cu(2)] are located adjacent to each 
other. The spins on the Cu(l) sites are in the nearly spin- 
singlet state, and a finite energy (dimer gap) is necessary 
to break the nearly spin-singlet state. The spins on the 



Cu(2) sites are almost polarized in the 1/2 plateau region, 
whereas they form AFLRO in low fields. Thus, unex- 
pectedly, coexistence of the nearly spin-singlet state and 
AFLRO appears, although the interaction between the 
Cu(l) and Cu(2) spins cannot be ignored. To our knowl- 
edge, Cu2CdB2 0g is the first material in quantum spin 
systems which shows such coexistence in an atomic scale. 
Future studies must address nuclear magnetic resonance 
(NMR) measurements and theoretical determination of 
a phase diagram in the spin system. 
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FIG. 1: (color), (a) Schematic drawing of Cu 2+ -ion posi- 
tions in CU2CCIB2O6. Red and blue circles indicate Cu(l) 
and Cu(2) sites, respectively. Red, green, blue, dashed 
black, and solid black bars represent Cu-Cu bonds 1, 2, 
3, 4, and 5, respectively, (b) An illustration of the spin 
system in Cu2CdB2 0e. The Hamiltonian is expressed as 

T~L = Ylil^^i,? ■ Si t 3 + J2(Si t l ■ Si, 2 + Si, 3 ' Si,i) + Js{Si,l • 

Si+1,1 + SiA ■ Si+1,4)]. Exchange interactions in the bonds 
1, 2, and 3 (Ji, J2, and J3, respectively) are AF and mainly 
determine magnetic properties. 



formed by the Cu(l) spins is small. As a result, the spin 
gap is expected to be small, because it vanishes when 
the effective interaction is zero. In the parameter regime 
corresponding to Cu2CdB2 06, the spin gap is too small 
to be detected in the scale of Figs. 2(a) and 3. Thus, it 
is expected that the spin gap is easily destroyed by weak 
3D couplings. 

[15] M. Oshikawa, M. Yamanaka, and I. Affleck, Phys. Rev. 
Lett. 78, 1984 (1997). 



FIG. 2: Magnetic susceptibility \{T) below 300 K in (a) and 
50 K in (b). The solid and dotted curves respectively indicate 
X (T) of Cu 2 CdB 2 6 and the QMC result of X (T) for the spin 
system in Fig. 1(b) with Ji = 160, J 2 = 38.8, and J 3 = 9.7 
K. The insets in (a) and (b) show respectively inverse x(T) 
and dx(T)/dT curves. The dashed line in the inset of (a) 
denotes the fitted line to inverse x(^) above 180 K. A single 
division of the vertical scale in the inset of (b) means 5 x 10~ 4 
emu/Cu mol K. (c) Specific heat C(T) of Cu 2 CdB 2 06 below 
50 K. The inset shows C(T) below 300 K. 



FIG. 3: Magnetization of Cu2CdB206 (solid curve) expressed 
as a spin value S = M(H) / q^bNa where /j,b and Na are the 
Bohr magneton and Avogadro's constant, respectively. Di- 
amonds, circles, and squares respectively indicate the QMC 
results of M(H) for total, Cu(l), and Cu(2) spins of the spin 
system in Fig. 1(b) with Ji = 160, J 2 = 38.8, and J 3 = 9.7 
K. Curves 1, 2, and 3 in the inset shows the QMC results of 
M(H) up to 140 T for total, Cu(l), and Cu(2) spins, respec- 
tively. 
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FIG. 4: (a) Schematic picture of the spin state in low fields. 
Ellipses suggest that the spins on Cu(f ) sites are in a nearly 
spin-singlet state with a dirtier gap. Solid arrows indicate that 
the Cu(2) spins form AFLRO, but the direction of the ordered 
spins has not been determined experimentally, (b) The spin 
state in the 1/2 magnetization plateau region. Solid arrows 
indicate that the spins on Cu(2) sites are almost polarized 
parallel to the applied field. 
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